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Abstract 



The next to leading order chiral corrections to the SU(2) x SU(2) Gell-Mann-Oakes-Renner 
(GMOR) relation are obtained using the pseudoscalar correlator to five-loop order in perturba- 
tive QCD, together with new finite energy sum rules (FESR) incorporating polynomial, Leg- 
endre type, integration kernels. The purpose of these kernels is to suppress hadronic con- 
tributions in the region where they are least known. This reduces considerably the system- 
atic uncertainties arising from the lack of direct experimental information on the hadronic 
resonance spectral function. Three different methods are used to compute the FESR con- 
tour integral in the complex energy (squared) s-plane, i.e. Fixed Order Perturbation Theory, 
Contour Improved Perturbation Theory, and a fixed renormalization scale scheme. We ob- 
tain for the corrections to the GMOR relation, 5 n , the value 5 n = (6.2 ± 1.6)%. This re- 
sult is substantially more accurate than previous determinations based on QCD sum rules; 
it is also more reliable as it is basically free of systematic uncertainties. It implies a light 
quark condensate (0|«u|0) ~ (0|dd|0) = {0\qq\0)\zGeV = (-267 ± 5 MeV) 3 . As a byprod- 
uct, the chiral perturbation theory (unphysical) low energy constant H 2 is predicted to be 
Hl{v x = M p ) = -(5.1 ± 1.8) x 10- 3 , or H r 2 {v x = Af„) = -(5.7 ± 2.0) x 10~ 3 . 
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1 Introduction 



The method of QCD sum rules [I] provides a successful analytical technique to determine hadronic 
parameters and form factors, as well as to extract from data fundamental parameters of QCD, 
such as quark masses and the strong coupling. The method is based on two fundamental pillars, 
(i) the Operator Product Expansion (OPE) of current correlators at short distances, extended 
beyond perturbation theory, and (ii) Cauchy's theorem in the complex energy (squared) plane, 
also known as quark-hadron duality. The extension of the OPE effectively accounts for quark- 
gluon confinement by parameterizing propagator corrections in terms of a set of vacuum matrix 
elements of the quark and gluon fields entering the QCD Lagrangian, i.e. the so called vacuum 
condensates. 

In this paper we concentrate on the pseudoscalar current correlator with pionic quantum num- 
bers, i.e. the two-point function involving the up- and down-quark degrees of freedom. Specifi- 
cally, we will study this correlator at zero momentum, which through a Ward identity leads to 
the SU{2) x 517(2) Gell-Mann-Oakes-Renner relation [2]- [5] 

if>s(P) = ~( m u + m d ) <0| uu + dd |0> = 2/| Ml (1 - 6 n ) , (1) 

where f n = 92.21 ± 0.14 MeV [S]-[Z]- Our goal is to determine the corrections to this relation, 
i.e. 5 n . In addition to its intrinsic importance to QCD, this quantity is of great interest to chiral 
perturbation theory (CHPT), as it is related to two of its low energy constants, as well as to 
the corrections to the SU(3) x SU{3) GMOR relation, i.e. at next to leading order [8j-[9j 

M 2 Ml 
5 w = 4-^(2L r s -H r 2 ) and 5 K = -^6 7T , (2) 

J TT 7T 

where Lg is a physical low-energy constant, and i7£ an unphysical one. The use of QCD sum 
rules to determine the chiral corrections to the GMOR relation has been fully justified in QCD 
in [5], |10j-[13] (for a recent reassessment see [H]). If one were to have direct experimental data 
on the hadronic spectral function (beyond the pion pole), the task of determining 5 n would be 
rather straightforward. Unfortunately, one only knows from experiment the mass and width of 
the first two radial excitations of the pion. This information is hardly enough to reconstruct the 
full hadronic spectral function. In fact, inelasticity, non-resonant background, and resonance 
interference are difficult, if not impossible to guess. Many attempts have been made in the past 
to model this spectral function, starting with threshold constraints from CHPT imposed on 
Breit-Wigner forms [15]. However, a serious systematic uncertainty is unavoidable. Historically, 
this has affected the accuracy of results from analyses based on the pseudoscalar correlator, 
such as e.g. QCD sum rule determinations of (light) quark masses [TJ, and estimates of the 
corrections to the GMOR relation [TT], [16]. Recently, a major breakthrough has been achieved 
by introducing analytic integration kernels in QCD Finite Energy Sum Rules (FESR) tuned to 
suppress the (unknown) contribution of the resonance sector [13], |17|-|18|. In this way it has 
been possible to reduce this contribution to a few percent of the total, which is made up by the 
pion (kaon) pole plus perturbative and non-perturbative QCD. Results for the up-, down-, and 
strange-quark masses, as well as the ratio of the strange to non-strange quark condensate are 
now available with basically no systematic uncertainty from the hadronic resonance sector. 

In [14] the correction 5^ was obtained indirectly by determining the ratio of the strange to the 
non-strange quark condensates, R su = (0|ss|0)/(0|uu|0), together with the CHPT estimate [TT] 
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L\{y x = M p ) = (0.88 ± 0.24) x 10" 3 . The integration kernel used in the QCD FESR was chosen 
as a second degree polynomial constrained to vanish at the peaks of the two radial excitations of 
the ground state pseudoscalar meson. The method of Fixed Order Perturbation Theory (FOPT) 
was used to compute the FESR contour integral in the complex energy (squared) plane. The 
results thus obtained are [Hj 

H r 2 {v x = M p ) = -(4.3 ±1.3) x 10" 3 and S n = (4 db 2) %. (3) 

In this work we determine 5^ directly from a QCD FESR involving Legendre-type polynomials 
as integration kernels. These kernels are tuned to reduce considerably the importance of the 
hadronic resonance contribution, thus reducing the systematic uncertainty arising from the lack 
of direct data in this region. This tuning is done in two ways, (i) as in |14j . |17]-|18]. i.e. using 
a second degree polynomial constrained to vanish at the peaks of the two radial excitations of 
the pion (local constraint), and (ii) using more general Legendre type polynomials with global, 
rather than local, constraints. This direct determination leads to a considerably more accurate 
result for i.e. an estimated total uncertainty of 25 %, rather than 50 % as in [T4], with both 
results being in agreement within errors. 



2 Finite energy sum rules and the Gell-Mann-Oakes-Renner re- 
lation 



We begin by introducing the light quark pseudoscalar correlator 

Ms = -q 2 ) = % J d^xe^ < 0| T(j 5 (x) j 5 (0)) |0 > , (4) 

where < 0| is the physical vacuum and the current density js(x) is 

j 5 (x) = {m d + m u )~d(x)i^u{x) , (5) 

and m u $ are the light quark masses. Invoking Cauchy's theorem in the complex energy (squared) 
s-plane, the correlator at zero momentum is given by 



1 f s ° Ac(s) If 
^ 5 (0)A 5 (0) = - / -^Im^s) ds+—i 

71 Js th s 2m f C n sol 



A5 S i/J5{s)ds, (6) 



where Sth = is the hadronic threshold, the integration path C( | sq | ) is a circle of (finite) 
radius so, and As(s) is an arbitrary, but analytic integration kernel. The radius sq determines 
the onset of the continuum, which is usually assumed to be given by perturbative QCD (PQCD). 
Hence, the first integral in Eq.(6) involves the hadronic spectral function, and the second integral 
the QCD correlator. Saturation of the hadronic integral with only the pion pole leads to the 
GMOR relation, Eq.(l) with S n = 0, after invoking the Ward identity 

^ 5 (0) = -[m u {^)+m d {^)]{Q\uu + dd\Q){^). (7) 

If quartic quark mass corrections to tp5(s) are neglected, quite legitimate for the up and down 
quarks, the condensate in Eq.(7) can be considered normal-ordered (for a discussion about these 
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corrections see [TU], [IS])- Both the quark masses and the normal ordered vacuum condensates 
depend on the renormalization point /x 2 , but their product does not. 



In the QCD sector the pseudoscalar correlator has been known for quite some time up to fourth 
loop order in PQCD |19j . and more recently the fifth order has been calculated |20j. Quartic 
quark mass corrections, as well as vacuum condensate terms are also known, but we found their 
contribution to be negligible in this application. In fact, this contribution is much smaller than 
the uncertainty due to the value of the strong coupling. Hence, the two-point function entering 
the contour integral in Eq.(6) can be generically written as 



/ \ i r / 2\ / 2\"|2 3 v > i (i) / 2\ ( C^sifJ^ 

^5{s)\pQCD = - [m u {n ) + m d {fi )J ps^^' I — - 



where 



i+i 

3=0 



In 



(i = 0, . . . 4) 



and the constant coefficients Cj,- are given in [19J-[20J. 



(8) 



0) 



The hadronic spectral function can be split into the pion pole term followed by the resonance 
contribution, i.e. 



1 



1 



-Im^{s) = 2tiM*8(s-Mt) + - Im ^{s)\ RES , 



7T 



TT 



(10) 



where the resonance contribution, Im ip5(s)\nES, involves the radial excitations of the pion, with 
the first two, 7r(1300) and 7r(1800), known from experiment [6]. This experimental information, 
though, is restricted to the values of the mass and the width. The resonance spectral function 
in this channel is completely unknown experimentally, so that its reconstruction is model de- 
pendent. This leads to a serious systematic uncertainty. In an attempt to reduce this model 
dependency, it was first proposed long ago [15] to normalize this spectral function at threshold 
using CHPT constraint. 
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TT 



9{s - 9M 2 



1 Mi 1 



*' 9 /2 2 8 ^ 4 
where the phase-space integral I\ps(s) is given by [2T] 



\PS 



(*) 



du \ 1 



I\ 



PS 



(11) 



5 + ^ 



(s - 3n + 3M^) 2 + 3 \(s,u,M*) 1 



u 



2 (s-M|)2 
+ 20 Mi 



+ 



(s - M*) 



3(u - Ml) -s + 9M% 



(12) 
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where 



X(l,u/s,M 2 /s) = 



u+M % y 



(VTi-M w ) 



In the chiral limit the phase space integral I\ps(s) reduces to the simple expression 

lim I\ps(s) = 3s , 

which leads to the well known Pagels & Zepeda result |22j 



- Im if>5(s)\, 

TT 



1 M 4 1 



3 /| 2^ ■ 

Using Eq.(ll) to normalize a model involving two Breit-Wigner forms leads to 



1 



TT 



Im ips(s)\ 



RES 



Im ^5(s)|, 



\BWi(s) + k BW 2 (s)} 

(l + K) 



where BWi(0) = 6^(0) = 1, with 



BWi(s) 



MfiMj +rf) 

(s _ M 2)2 +M 2 r 2 



(t = l,2) 



(13) 
(14) 

(15) 
(16) 
(17) 
(18) 



and ft is a free parameter controlling the relative weight of the resonances. Imposing that the 
area under the spectral function in zero width be equal to the area under the spectral function 
in finite width, in the limit r — >■ 0, will then determine the value of the leptonic decay constant 
of the first radial excitation, f ni , if only the first resonance is used. In fact, using the relation 



lim 



r^o (s - Ml) 2 + Mfrf Mil?! 
gives for f m the prediction 



Ml 1 



Mi 



it 5(s - Ml) , 



1/2 



(19) 



(20) 



IW1 f v 16 V67r 3 r\ 

Notice that the correct chiral behaviour of f ni is obtained, i.e. f Wl = O(Ml), as / 7ri must vanish 
like M 2 in the chiral limit because tt\ is not a Goldstone boson (the chiral symmetry remains 
SU(2) x SU (2) in the presence of radial excitations of the pion). Using in Eq.(20) Mi ~ 1.3 GeV, 
and ri ~ 200 - 400 gives 



2 MeV . 



(21) 



In order to reduce the model dependency of the hadronic resonance parametrization it was 
proposed recently [H], [Ej-[l8] to choose an integration kernel in the FESR, e.g. in Eq.(6), such 
that it vanishes at the peak of each resonance 

A 5 (s) = 1 - a s - ax s 2 , (22) 

so that A 5 (Ml) = A 5 (M|) = 0, which fixes a = 0.894 GeV 2 , and a x = -0.179 GeV 4 . This 
simple kernel does achieve a substantial reduction of the systematic uncertainty arising in the 
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hadronic sector. In fact, the resonance contribution to the up-, down- and strange-quark masses 
is at the level of a couple of percent, well below the error due to the uncertainty in the strong 
coupling. The optimal choice of integration kernel is most likely application dependent. For 
instance, pinched kernels in FESR for the vector and axial-vector correlators [23J improve con- 
siderably the agreement between data and the first two Weinberg sum rules. However, a com- 
bination of Eq.(22) and pinched kernels does not work as well as Eq.(22) on its own |18j . 



A generalization of Eq.(22) in the form of a n-th degree Legendre type polynomial has been 
used in the past to determine meson coupling constants [24], heavy quark masses |25j . and to 
extract chiral vacuum condensates from data, as well as to determine the counter term of the 
0(p 4 ) CHPT Lagrangian [26J. We shall also use such a kernel here in order to compare with 
the simple form, Eq.(22). We introduce the general polynomial 

n 

P n (s) = Y, a m s m , (23) 

m=0 

where the coefficients are fixed by (i) imposing a normalization condition at threshold 

P n {s = Ml) = 1 , (24) 
and (ii) by requiring that the polynomial P n (s) minimizes the contribution of the continuum in 



the range 



(fc+i) 
'th > s o 



in a least square sense, i.e., 

so 

s m P n (s) ds = for m = 0, . . . n - 1 . (25) 

th 

The meaning of the threshold s^ 1 ^ * s as foh° ws - The hadronic resonance spectral function 
can be split into the sum of a term involving the known pionic radial excitations, and a term 
involving unknown resonances, non-resonant background, etc., i.e. 



Im^f \ s ) = J^ImiP^is) + 6(s - s[* +1) ) Im4>f +1 \s) , (26) 

i=0 



where the upper index k stands for the number of resonances explicitly included in the sum. The 
polynomials obtained in this way are closely related to the Legendre polynomials, so that Eq.(25) 
is exact on account of orthogonality. Specifically, denoting by V n {x) the Legendre polynomials, 
their relation to the P n (s) above is 



where the variable x(s) is 



( ) _ Vn{x{s)) 



defined in the range x G [—1,1] for s E [s{:^ +1 \so]. The normalization of the V n (x) is such 
that Vi(x) = x, Vi{x) = (5x 3 - 3x)/2 , etc. 
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3 Integration in the complex s-plane 



The FESR determining ips(0), Eq.(6), can be written as 

^ 5 (0) A 5 (0) = 2 ft Ml A 5 (M2) + S 5 (s )\ RES + 8 5 (s )\ QCD , (29) 



where 



and 



1 f s ° A 5 (s) 

$5(so)\res = - Imip 5 (s)\ RE s ds , (30) 

hul s 



h(so)\QCD = tt~- / A5 S ^{s)\qcd ds , (31) 
2m fc(\s \) s 

and equivalent expressions for the Legendre type polynomial kernels, where As(s) becomes P n {s) 
as defined in Eqs. (23)-(28). The contour integral in the complex s-plane can be evaluated in 
the framework of Fixed Order Perturbation Theory (FOPT), or Contour Improved Perturbation 
Theory (CIPT)[27j. In FOPT the coupling and the quark masses at a scale sq are considered 
fixed (constant) so that only logarithmic terms contribute to the integral. The renormalization 
group (RG) summation of leading logs is only carried out after the contour integration by setting 
{j? = —sq. In this case the integration can be done analytically. In fact, defining 

J(q,k) = ^-J s«(ln^f) k ds, (32) 

where k = 0, 1, ...5 and q = 0, 1, .., n, one can derive the master formula 

p=0 1=0 1 vi < j \ r / 

Explicit expressions for 5^{sq)\qcd may be found in [13]. In the case of CIPT the strong 
coupling and the quark masses are running and the RG is implemented before integrating. The 
FESR then must involve the second derivative of the pseudoscalar correlator which is introduced 
through the identity 



dsg(s)Ms)=f ds [F(s) - F(s )} ^'(s) , (34) 
C(\s \) JC(\so\) 



where 

F(s) = I ds' 



ds g(s ) - / ds"g{s") 



(35) 



and g(s) is an arbitrary function. This identity is easily proved by integrating by parts the right 
hand side to obtain the left hand side. We choose g(s) = A5 (s)/s, with As(s) given in Eq.(22). 
The function F(s) then becomes 

F(s) = s ln{s/s ) - (l + a s + ^ a x s 2 J s + ^ s 2 + ^ s 3 . (36) 
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The running coupling and the running quark masses can be computed at each step in the 
integration contour by solving numerically the corresponding RG equations (RGE)[28j-[29j. The 
running coupling satisfies the equation 

da s (—s) _ ^ _ ^ fl _ n f 



p{a 8 ) = ~^2Pn a s (-s) N+2 , (37) 



ds 

N=o 

where a s = a s /ir, and for three quark flavours /3o = 9/4, fi\ = 4, fa = 3863/384, = 
(421797/54 + 3560C(3))/256. In the complex s-plane s = soe lx with the angle x defined in the 
interval x S (— 7r,7r). The RGE then becomes 

^M = - l J2PNa s (x) N ^. (38) 

X N=0 

This RGE can be solved numerically using e.g. a modified Euler method, providing as input 




3.8 4.0 4.2 
So (GeV2) 



5.0 



Figure 1: Results for ^ 6 (0) in units of 10~ 4 GeV 4 function of so in FOPT and neglecting the 
hadronic resonance contribution. Curve (a) corresponds to a s (M 2 ) = 0.335 (A = 365 MeV), 
and curve (b) to a s (M?) = 0.353 (A = 397 MeV). 



a s (x = 0) = a s (— sq). Next, the RGE for the quark mass is given by 



s dm(—s) 
m ds 



M=0 



(39) 



where for three quark flavours 70 = 1, 71 = 182/48, 72 = [8885/9 - 160C(3)]/64, 73 = 
[2977517/162 - 148720 C(3)/27 + 2160 £(4) - 8000C(5)/3]/256. With the aid of Eqs. (37)-(38) 
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the above equation can be converted into a differential equation for m[x) and integrated, with 
the result 

m(x) = m(0) expj -i dx' 7a/ [a s (x)] M+1 } , (40) 

"'° A/=0 

where the integration constant m(0) is identified as the overall multiplicative quark mass in the 
expression for the pseudoscalar correlator, i.e. m = [m u (so) + md(so)]. 



> 
o 




3.8 4.0 4.2 
So (GeV2) 



5.0 



Figure 2: Results for 1/15(0) in units of 10 4 GeV 4 as a function of sq in CIPT and using a two- 
resonance parametrization of the hadronic spectral function normalized at threshold according to 
CHPT. Curve (a) corresponds to a s (M?) = 0.335 (A = 365 MeV), and curve (b) to a s (M*) = 
0.353 (A = 397 MeV). 



4 Results 



We begin by discussing the results of the first method, i.e. using the integration kernel Eq.(22). 
This choice reduces considerably the resonance contribution to ^5(0), which becomes smaller 
than the uncertainty due to the error in the strong coupling, both in FOPT as well as in CIPT. 
In fact, the difference between including and not including the resonance sector is at the level of 
less than 1% in the result for ^5(0). We use the recent high precision determination [30] of the 
strong coupling a s (M^) = 0.344 ± 0.009, which corresponds to a QCD scale in the MS scheme of 
A = 365 — 397 MeV. This value is in perfect agreement with the result of a more recent analysis 
using all present experimental and theoretical knowledge [31], i.e a s (M^) = 0.342 ± 0.012; the 
tiny difference between these two values of the coupling have no impact in our final results. 
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For the quark masses we use the most accurate determination using the kernel Eq.(22) [18] . 
essentially free of systematic uncertainties, m^(2GeV) = 2.9 ± 0.2 MeV, and m^(2 GeV) = 
5.3 ± 0.4 MeV. The uncertainties are mostly due to the strong coupling. These running masses 
correspond to the following invariant masses (the integration constants in the RGE) m u = 3.8 — 
3.9 MeV and = 6.4 — 7.2 MeV. The results for 1/15(0) in FOPT, and neglecting completely 
the hadronic resonance contribution is shown in Fig.l for a s (M 2 ) = 0.335 (A = 365 MeV) 
(curve(a)), and a s (M 2 ) = 0.353 (A = 397 MeV) (curve(b)). These results are fairly stable in 
the wide region s = 3.0 - 5.0 GeV 2 , giving ^ 5 (0) = 3.12 - 3.14 x 10~ 4 GeV 4 for A = 365 MeV, 
and i/i 5 (0) = 3.08-3.09 x 10" 4 GeV 4 for A = 397 MeV. Taking now the resonance contribution 
into account, and using a two-resonance parametrization normalized to CHPT at threshold, as 
described in Section 2, 1/15(0) increases by less than 0.6% for each value of A. The principal 
source of error is then that from a s , from the sq dependence, and obviously from the unknown 
six- loop order in PQCD. Combining these values of ^5(0) leads to the following result for S n 

SttIfopt = (6.5 ± 0.9)%. (41) 

In Fig. 2 we show 1/15(0) obtained in CIPT with the two-resonance parametrization of Section 2 



3.18 - 
3.17 - 
3.16 




So (GeV2) 



Figure 3: Results for 1/15(0) in units of 10~ 4 GeV 4 as a function of sq in a fixed fi 2 = 4 GeV 2 
scheme, and for the central value a s (M 2 ) = 0.344. Curves (a), (b), (c) correspond to polynomial 
degrees n = 4, 5, and 6, respectively. 



for both values of a s . Once again due to the presence of the integration kernel, the systematic 
uncertainty due to the resonance contribution is much smaller than the uncertainty due to the 
value of a s (M 2 ) and the s dependence. The result is ^ 5 (0) = 3.07 - 3.10 x 10" 4 GeV 4 , 
and ip 5 (0) = 3.12 - 3.15 x 10" 4 GeV 4 , corresponding, respectively, to a s (M 2 ) = 0.353 (A = 
397 MeV), and a s (M 2 ) = 0.335 (A = 365 MeV). These results translate into 

S^lciPT = (7.0 ± 0.8)%, (42) 
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in very good agreement with the FOPT result, Eq.(41). 

Next, we use the Legendre-type integration kernel, Eq.(23), together with the global constraint, 
Eq.(25). The pseudoscalar correlator is known to be affected by the rather slow convergence of 
its PQCD expansion. This is due to the monotonically increasing coefficients multiplying the 
higher powers of a s . As a consequence of this, the convergence of the results will depend crucially 
on the renormalization prescription being used. We favour a pragmatic a posteriori approach. 
The preferred scheme should be the one that yields the most stable results with respect to 
two parameters, the degree of the polynomial P n (s) and the duality radius so- As far as the 
renormalization scheme is concerned there are two obvious choices related to the renormalization 
point in the asymptotic expansion: we can either fix /x 2 = s (CIPT) to avoid large logarithmic 
terms in the perturbative expansion or, since the logarithms entering our analysis are not really 
large yet, use a fixed reference scale, say /x = 2 GeV. In the case of the Legendre-type kernels 
we find that a fixed value of /i provides better stability. In Fig. 3 we show the results for "05(0) 
corresponding to n = 4,5,6, and the central value a s (M 2 ) = 0.344. Using values of a s in the 
range a s (M 2 ) = 0.335 — 0.353 introduces an uncertainty in 5 n of 5%, and changing the fixed 
scale /d 2 in the range /x 2 = 2 — 50 GeV 2 produces a larger 18% uncertainty, leading to 

6^ = (5.6 ± 1.1)% , (43) 

compatible with the FOPT value, Eq.(41), as well as the CIPT result, Eq.(42). The dominant 
uncertainty is that due the /x dependence, as indicated above. If one were to use somewhat 
lower values of the strong coupling, e.g. a s (M 2 ) = 0.322 ± 0.02 in FOPT [32], the result would 
fall within the error band in Eq.(43). Combining all three results conservatively by using the 
whole spread of values, rather than averaging and adding errors in quadrature, leads to the final 
combined result 

5 n = (6.2 ± 1.6) % , (44) 

with an uncertainty about 50 % smaller than that in Eq.(3) [13] • This value of 5 n is somewhat 
higher than earlier results based on PQCD to NLO (21]: 5 n = (4.0 ± 1.0)%, and to NNLO [33]: 
5 n = (3.5 ± 1.0)%, together with model dependent hadronic spectral functions (thus having a 
far greater impact on the result than in the present determination). Using the result, Eq.(44), 
in Eq.(l), together with the recent most accurate quark mass value [TBj \{jn u + m c i)\2GeV = 
4.1 ± 0.2 MeV, and assuming (0|mt|0) ~ (0|dd|0) = (0|gg|0), we find 

<0|gg|0)| 2G eV = (-267±5MeV) 3 . (45) 

The result, Eq.(44), implies, using Eq.(2) together with [UJ L r 8 (u x = M p ) = (0.88±0.24) x 10~ 3 , 
the following prediction 

H r 2 {v x = M p ) = -(5.1 ± 1.8) x 10~ 3 , (46) 

in reasonable agreement within errors with the estimate |11] H^^x = M p ) = — (3.4±1.5) x 10~ 3 , 
as well as the determination in Eq.(3) [14] . Using, instead, the lattice QCD result [33] L%(y x = 
M v ) = (0.58 ± 0.09) x 10" 3 gives 

H r 2 (v x = M v ) = -(5.7 ± 2.0) x 10" 3 . (47) 

In closing, a few remarks should be made concerning the input used for the strong coupling. 
First and foremost, the issue addressed in this paper is related to physics at scales closer to 



10 



the r-lepton mass than to e.g. the Z-boson mass. For this reason, we have chosen to use 
FOPT and CIPT direct determinations of a s (M T ). The three sources used |30j-|32j lead to 
values of 5 n within the error band in Eq.(44). Additionally, determinations of a s (M T ) are based 
on a transparent, model independent procedure which uses as input experimental data on r- 
decays. This should be contrasted with determinations of the strong coupling at the scale of 
the Z-boson mass, which involve far more complicated dynamics. The latest determination of 
a s (Mz) |35| gives results which if extrapolated down to the r-mass would somewhat disagree 
(at the one standard deviation level) with the direct determinations of a s (M T ). In fact, the 
authors of [35] find a s (M z ) = 0.1135 ±0.0010, implying a s {M T ) = 0.305 ±0.008 to be compared 
with the range a s (M T ) = 0.320 — 0.353 from direct determinations from r-decay data in FOPT 
and CIPT. A previous world average at the Z-boson scale [36] gives the slightly higher value 
a s {M z ) = 0.1184±0.0007, or a s (M T ) = 0.3186±0.0058. The reconciliation between these results 
and those at the r-scale is beyond the scope of this work. Assuming that the problem does not 
lie with the r-data, our final result, Eq.(44), incorporates the current precise information on 
a s {M T ) from both FOPT and CIPT. 



5 Conclusions 



In this paper we have performed a direct determination of the light pseudoscalar correlator at 
zero momentum, Vs(0)> using QCD FESR involving polynomial integration kernels tuned to 
suppress the experimentally unknown hadronic resonance spectral function. This suppression 
is successfully achieved, as in the determination of the strange quark condensate [2], and of 
the light quark masses |17j-|18|. In fact, the uncertainty due to the resonance contribution to 
the FESR is much smaller than the uncertainty due to a s , and the so dependence. Two types 
of integration kernels have been used in the FESR. First, a simple second degree polynomial 
vanishing at the peaks of the first two radial excitations of the pion (local constraint). Second, 
Legendre type polynomials of varying degrees with global constraints. The fact that results from 
both methods are consistent with each other provides strong support for this procedure. The 
integration in the complex s-plane has been performed according to three methods, i.e. FOPT, 
CIPT, and a fixed renormalization scale scheme. Results from the different integration kernels, 
and the different integration procedures are consistent with each other and lead to a substantial 
reduction of the uncertainty in the corrections to the GMOR relation. This uncertainty is mostly 
due to the value of the strong coupling, the so dependence, and to a lesser extent due to the 
integration procedure and the nature of the integration kernels. The systematic uncertainty 
from the resonance sector is basically eliminated by these kernels. As a byproduct, the value of 
the light quark condensate, and the CHPT low energy constant have been obtained. 
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